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We discuss a large class of phenomenological models incorporating quantum gravity motivated 
corrections to electrodynamics. The framework is that of electrodynamics in a birefringent and 
dispersive medium with non-local constitutive relations, which are considered up to second order in 
the inverse of the energy characterizing the quantum gravity scale. The energy-momentum tensor, 
Green functions and frequency dependent refraction indices are obtained, leading to departures from 
standard physics. The effective character of the theory is also emphasized by introducing a frequency 
cutoff SI. The analysis of its effects upon the standard notion of causality is performed, showing 
that in the radiation regime {Q,R » 1) the expected corrections of the order {u/Q.)" get further 

I suppressed by highly oscillating terms proportional to smiyiR), cos{QR), thus forbiding causality 

, violations to show up in the corresponding observational effects. 
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. I. INTRODUCTION 

' Lorentz covariance in local inertial frames is a well established symmetry at the energies of present day experiments. 

(N 

However, its validity at high energies is subject to test. Possible Lorentz invariance violations may arise from dynamical 
<N . . . 

^ modifications induced by quantum gravity (QG). The effects of such violations in the range well below the Planck 

m: ... finnnn 



This issue is closely linked 
61], concerning Lorentz and CPT 



\^ . energy {Ep 10 GeV) have been recently the object of intense scrutiny 

o . 

0^ ■ to theoretical and experimental research, based on the Standard Model Extension 
\^ violations [7] . Heuristic loop QG derivations of such effects [2|, |3| make it clear that a better understanding of the 

o 

■ corresponding semiclassical limit is required 8]. String theory has also provided models for explaining such QG 
induced corrections [9|. Moreover, effective field theory models have been constructed that include higher dimension 

^ ... n . n 

Lorentz invariance violating (LIV ) operators [lO|. Synchrotron radiation arising from the model in [10] has been 
. !_| \ extensively analyzed in Ref . [ll| . These effective theories use a reduced number of degrees of freedom to describe 
^ ■ the physics at a low energy scale, ignoring the detailed dynamics inherent to Planck energies. In other words, if QG 
dominates at a scale Eqq, usually considered of the order of Ep, a corresponding low energy effective theory can be 
visualized as an expansion in powers of ^ Eq^, truncated at a given finite order. In this way it will be a good 
description, hopefully simpler than the original one, for energies E <C Eqq. This restricted validity relaxes some of 
the constraints usually required for physical theories, such as renormalizability. Stability and causality, perhaps of 
more essential status than the Lorentz symmetry itself, are assumed to remain valid at the low energy regime [l^ . 
Nevertheless, fine tuning problems arise when considering radiative corrections d2], which can be circumvented by 
extending the notion of dimensional regularization |l^. 

In fact, one of the possible manifestations of QG at low energy is the appearance of correction terms related to the 
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scale Eqq in the standard particle prop agation and interaction properties. The most direct interpretation of such 
corrections, though not the only one ^15], is in terms of a spontaneous breaking of Lorcntz covariancc at high energies. 
If this is so, the effective theory will be covariant under Lorentz transformation between inertial frames (passive or 
observer transformations), and the observable Lorcntz symmetry violations will be associated with rotations and boosts 
of the fields in a given inertial frame (active or particle transformations). In this case the space-time coordinates at low 
energy remain commutative. We focuse here on a general analysis of QG effects in electrodynamics in these models, 
where we can introduce the full usual mathematical framework of field theory, specially Fourier transformations. QG 
may also modify the space-time itself so that the coordinates become noncommutative, as in the case of Double Special 
Relativity models, for example. In this situation, which will not be discussed here, ordering ambiguities preclude a 
direct use of such transformations. 

Considering now the electromagnetic field, the models proposed to describe low energy effects of QG can usually be 
expressed in terms of modified dispersion relations, with a polynomial dependence in energy and momentum. Such 
modifications include standard Lorentz invariance violations as well as possible extensions of Lorentz covariance (l7| . 
Most of these approaches can be unified in the description 

ik-D = 47rp, k-B = 0, (1) 

kxE-wB = 0, ik X H + icjD = 47rj, (2) 

where the auxiliary fields 

D' = a'^Ej + p'^Bj , W = P'^Bj + a'^ Ej , (3) 

are such that the coefficients a''^ , , p*-' and ct*-' depend on the energy cu and the momentum k of the electromagnetic 
field. These equations correspond to a higher order linear dynamics. Equations ([1]) and ([2|) strongly resemble the usual 
description for an electromagnetic field in a medium [isl , where the fields D and H are characterized by constitutive 
relations of the form ([3]). In terms of electrodynamics in media, we can interpret the low energy QG corrections in 
terms of a dispersive bianisotropic media. From a heuristic point of view, as shown below, these effective media are 
non-local in space and time, which can be interpreted as a footprint of the granularity induced by QG. 

Strictly speaking, the effective models are characterized by Eqs.((l][2]), but under certain restrictions it is also possible 
to pose an action from which they derive. Although not essential, this is a useful approach to visualize general features 
of the dynamics. Let us recall the Lagrangian for the electromagnetic field in a local medium 

L = -iF^.x'^''""^'^^a/3 - 47^J^^^ (4) 

where F^^, = d^Ai, — di^A^, and contains the information about the medium. This structure warrants gauge 

invariance and hence charge conservation. The dynamics is given by the equations of motion 

d^H'^-' = Anf , (5) 



together with the constitutive relations 

^M<> ^ ^M["/3]^^^, (6) 

Defining the electric and magnetic fields as pQi = Ei and Fij = —CijkBk respectively, and the corresponding compo- 
nents oi H'^'^, = and iJ'-* — — e'^'^iJ'^, the constitutive relations become 

In our notation the components of any three-dimensional vector V are given by those with subindices Vi. 

The equations of motion incorporating QG corrections acquire a similar form, but with one important difference 
arising from the nonlocal character of the effective medium. The ;^[^''l[a/3] jg now a non-local tensor, such that 

instead of ([4]). As usual if F^^, and L are real, the reality of ^'^'^i'"''! (a;'^ — x'^) and subsequently of 

m%x)^ J d^ix^^^'^^^^^x" - i") F^ffiS:"). (9) 
are implied. Writing ^['^'''["'^'(a;'^ — 5°') in terms of its Fourier transform 

we can easily demonstrate that ([8]) can also be written a 

L = -^F^A^n [x'^''""''l(*a.) F^^ixn] , (11) 

with being a derivative operator. In terms of the Fourier transform, the reality of ^['''''["'^l (x'^ — x'^) is stated 

as 

(12) 

which holds similarly for the transformed fields F^pik'^) and Hapik'^)- If x''^''""''' is symmetric, in the sense that for 
each set of index values (/i,i^,a,/3) (no sum with respect to repeated indices) 

d^x F^. (x['"'""''1f„/3) = 1 d\ F^p (x'"^"'"'l^^p.) (13) 

is satisfied, then it is possible to perform integrations by parts making the equations of motion of the same form 
as in the usual non-operator case. Thus the Fourier transform of the equations of motion and constitutive relations 
acquire the structure of ((TJ21) and ^ respectively. In the following we assume that this property is satisfied. When 
the components of x''^*''!"''' do not correspond to a standard Lorentz tensor, this Lagrangian describes a model where 
the Lorentz symmetry is broken by the medium. 



We use this approach, where the QG modifications are described phenomenologically by constitutive relations, to 
discuss the main properties of QG induced effects in electrodynamics. A low energy expansion is developed in terms 
of the parameter ^ ~ ^qg- Working to order allows us to present jn the form 

where the constant coefficients ^ ■^^''^l^^'^'^I"^! are proportional to respectively. They are antisymmetric 

in the indices inside square brackets and symmetric in the indices inside curly brackets. In this way we are considering 
a Lagrangian depending up to third derivatives in the basic electromagnetic potential A^. Moreover, the integration 
conditions (|13p require the following symmetry properties 

[p.u]lal3] _ [al3][t^u] [t,u]e[al3]_ [a/5]0[MH MW}[a/3] _ [q/3]W}M (.r^ 

AO — Xo ' Al — ~ Al J A2 — X2 ■ l^Oy 

Once the coefficients of the constitutive relations have been promoted to derivative operators we obtain the relations 

^[o.][o,] ^ 1^.,^ ^ (16) 

X'""""-''' = +lemmij'^, x'"^"""' = +^efci»/9*^'e,mn, (17) 

by comparing Eqs. ([7]) and The expansion induces the corresponding form in the coefficients of the 

constitutive relations 

a-' —aQ + t.ai + 1. 0^2 o^oe, (tSj 

and similarly for p*^, ct*-' and in an obvious notation. Here ag'"''*, Q;^*"'^'^and a'^'''^^^ are constant coefficients. 

To analyze the propagation of the fields and to define the corresponding refraction index, the dependence of the 
constitutive relations on lu and k has to be made explicit. To achieve this we first expand the coefficients of the 
constitutive relations in space derivatives, maintaining covariance under rotations. Considering that these models can 
be understood as perturbative descriptions in terms of the parameter ^ we have, up to order 

a'' = aoidtW' + ai{dt)ie''''dr + «2(9t)l'5'5^ (19) 

with analogous expansions for , ct*-' and p*-' . Here ua, Pa, o^a, and pA, with A = 0,1, 2, are 5*0(3) scalar operators. 
This approach can be generalized to models with preferred spatial directions. The symmetry of x in Eq. (jl3[) implies 
that the terms in a*-', p'^and /3'^ with an even number of derivatives are symmetric under i ^ j, while the terms with 
an odd number of derivatives are antisymmetric. In the case of p'^ and a^^ Eq. (|13p leads to 

PA = -CTA- (20) 

Furthermore, we can also consistently expand the coefficients a a, Pa, o'a, and pA in powers of 1^ , according to 

Ca ^ Cao + CAiidt + CA2^^9f , (a = {aA, Pa, Pa, <Ja}, (21) 



with (ao, Cai, Ca2 being constant coefficients. This expansion leads to the identification of the following coefficients 

"2 -'7 ^00, P2 — V P02^ P2 — V P02, "2 - P2 " P2 —LI: 

^fe)mn ^ 1 ^^.,m^j„ _^ ^jnjjm^ ^^^^ /J^^^'"" = 1 + §'"5^"") (320, pf^''"' = ^ + ^m^jm) ^^O- (22) 

The above partition is consistent with the requirement of covariance under rotations. We have taken apo = /3oo = 1 
and Poo = cqo = to recover the usual vacuum as the background for ^ = 0. 

The fact that this theory does not hold at high energies will be coded by cutoff 57 ^ Eqg- We provide a general 
description of such modified electrodynamics including expressions for the equations of motion, the energy-momentum 
tensor and the Green functions as well as the corresponding refraction indexes, up to second order in ^. 

II. EQUATIONS OF MOTION 

Equations p9p . together with the corresponding ones for the remaining coefficients of the constitutive relations give 

D = (ao + aaC^k^^ E - (^cto + iai^u^ B + [icjii + wa2l^) (k x B) , (23) 
li={Po- iujcji^ B - i/3i| (k X B) + ctqE. (24) 

In the approximation to order here considered, we have ^^fc^ c± ^^w^ and we can write 

D = di(cj)E + id2(^^)B + d3(w)l(k X B), (25) 
H = /ii(w)B + i/i2(w)E + ih^{uj)i (k X B) , (26) 

where the functions (ii(tj) and hi{uj) depend only on oj and admit a series expansion in powers of ^w, characterizing 
each specific model. From Eqs. (IT][2]) we get the equations for E and B 

idi(k-E) =47rp(tj,k), (27) 
diE + (/igfc' - g{Lo)) |B + (wd3^~+/ii) {ik X B) = 47rj (a;, k) , (28) 



where we denote 

{d2 + h2)uj^ g(uj)l (29) 

The expressions p3ll24p indeed indicate that the above combination is of order ^. We thus see that in fact there are 
only three independent functions of uj and k which determine the dynamics 

P = di , Q = hi+ ujdii, R = {hsk^ - ff (^)) I (30) 



Using the homogeneous equation = k x E that yields w (k x B) = (k • E) k— fc^E, and charge conservation 
cjp — k • J = 0, we decouple the equations for fields E and B. Finally, we introduce the standard potentials <i> and A 

B = ik X A, E = iojA-i'k<^. (31) 

in the radiation gauge, k • A =0, in which case we have 

<^ = 4:n{k^Py^ p, (32) 

(fc^Q - uj^P) A + ii? (k X A) = 47r j - (^j ■ k^ k = inir, (33) 

from Eqs. (j27ll28p . The presence of birefringence depends on the parity violating term proportional to R. It is clear 
that a diagonalization is obtained in a circular polarization basis. Decomposing the vector potential and the current 
in such a basis 

A = A++A-, jT=jJ+jT> (34) 
and recalling the basic properties k x A+ = — «A+, k x A^ = iA~ , we separate (1551) into the decoupled equations 



[k^Q-uj^P + XkR]A^=ATrj^, X = ±l. (35) 
In terms of the basic functions introduced in the constitutive relations (|25ll26p , the factor in ([55]) is rewritten as 

fc^Q - Lu'^P + Afc i? = A/i3/c^|+ (hi + dsLU^) - Xgik ~ diuj"^ . (36) 



This is the key expression to obtain the Green functions and the refraction indices. 

III. GENERALIZED ENERGY-MOMENTUM TENSOR 

Any application of this modified electrodynamics related to radiation and its properties requires the construction 
of the corresponding energy momentum tensor. This section is devoted to such a construction. The theories under 
consideration are of higher order in the field derivatives, and thus call for an extension of the standard Noether 
theorem. The manipulations are highly simplified by proceeding in a covariant notation, the point being that the 
tensorial operator ^^'^'^l'"'') is constructed in a given reference frame and satisfies only passive Lorentz covariance. We 
assume that active Lorentz invariance is violated while active translation invariance is maintained, so that there is an 
energy momentum tensor given by the Noether theorem. 

Before constructing this tensor in our particular case we recall the general formalism for a Lagrangian including up 
to three derivatives in the fields. This can be useful also when consistently including dimension five and six operators 
in the matter field coupled to the above electrodynamics. We start from an action of the form 



(37) 



where we consider ^a, ^a^i, d,y^A,fj., di^dp^A.^L, to be independent fields, i.e. at this level we take for example that 
^v^A,^l 7^ dij,^A,u- Applying the standard action principle one derives the Euler-Lagrange (EL) equations 



Assuming that translations generated by x'^ = a;'' + are a symmetry of the action ([57)) . Noether's theorem leads 
to the energy momentum tensor 

= -51 L + [J^-dJ ,,J^] + d,d. ^ 



d^A,r ydid,^A.r)J ' \d {dpO,^ A.r) , 



whose conservation drT'^^ = can be directly verified via the equations of motion 

Next we apply the above general results to our Lagrangian pT|) together with the realization (HU, where $^ = A^. 
The corresponding derivatives are 

- U, jj — - -i'^L,^Xo - oer^fs - -Xi daOet^p, (4U) 

The equations of motion outside the sources can be written as 

= (42) 

where 

-{i^)- inih) - (ro) - 

Now let us consider the energy-momentum tensor p9p . Let us observe that the second term in this equation is 
precisely — ^a,o-^^^" which is not directly gauge invariant. It can be rewritten as 

- = -F„^H^^ - A„^^H^^ = -F„^H^" - 9„ (^,F™) , (44) 

by using the equations of motion. The last term is identically conserved and does not contribute to the corresponding 
charges. The remaining contributions are 

d{drA^,,)~^''[d{dpdrA^,,))^~2''' ^^^+2^2 dpFo^^ 

which are naturally gauge invariant. The final gauge invariant, non-symmetric energy- momentum tensor is 
T; = -5lL- F,,,H^^ - \ (a.F^„) xT^^^^'^^Fe^ 
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A direct but rather long calculation allows to verify the conservation SrTJ" = 0, via the equations of motion (|42p . 
To express the energy-momentum components in terms of the fields E and B we use the following (3 + 1) splitting 

^l-ivXl ^aP — 'iWoiXl ^Om + ^Xl [yyos'-Jmn — yVmn(-^Os\ + t^ijXl '-'mm, 

T/T/ ^jA'H{-rp}[a/3]7-r /I U/ ^,lOi]{Tp}lQm]^^ , 9, JOs]{Tp}[mri] j-^ 1 U/ TT 1 1 M/ ,,1^ 1 f'^P} ['"H r r (Ar,\ 

^^ii'X2 Uaf3 — ^yVoiX2 Uam + ^X2 [W QgU mn + VV mnU as\ + VV ijX2 Umm, I ) 

for antisymmetric fields W^i^, Ua/s- We also recall the relation 

-L~ Foo.H''" = i (E • D + B • H) , (48) 

which is useful in calculating the energy density. 

Let us illustrate the above construction by writing the energy density u ~ Tq and the Poynting vector Si — Tq in 
terms of the fields E and B to first order in ^: 



S = E X B+i^aioE x ^tE + |^io 



i {dtB) xB - Ex (V X B) 



~a^ol\dt [E X B] , (49) 



u = ^ (E2 + B2) - i/Jio^B • V X B - iaiolE • V X E - iaio^V ■ [E x B] . (50) 
The terms proportional to ctio correspond to the liberty of modifying the energy momentum tensor as 

f:^T:+dpVi^'>\ i/M^-^jp-], (51) 

which was previously used in Eq. pi)) . In the (3 + 1) partition the above means 

u = u - V ■ Q, S = S + ^iQ + V X W, (52) 

with Qi — vj'"' and Wi — |ey7c^'"''^' • The last terms in Eqs. and ([50)1 correspond to the choice 

Q = iaioaExB), W = 0. (53) 

Thus, the contributions proportional to ctio can be eliminated and one recovers the corresponding expressions that 
can be obtained directly from Maxwell's equations. Using the fields E , B, D and H, together with the equations of 
motion in vacuum and the freedom given by the energy-momentum tensor transformation (|51|) , these magnitudes 
can also be written to first order in <^ in a much more compact form as 

S-i(ExB + DxH)-i (aio - /3io) ^Bx (V x E) - Ex (V x B)) , (54) 
M= i(E-D + B-H). (55) 



IV. GREEN FUNCTIONS 



sis is 



The exact retarded Green function for the potential A in the circular polarization basis 



+ i'>^<^irk-r , (56) 



where R = r — r', ki = A:i/|k|, k = |k|, and G (w,k) is obtained from Eq. (l35l) . 

Taking the analytic continuation uj —^ lu + ie to obtain the causal Green functions, only the poles in the upper half 
plane of k make a contribution to the integration. By successive rescalings, the denominator in Eq. (|57|) can be 
written in a more convenient form 

Qfc2 - _^ ^i^j^ ^ -nluj^QaiM^ - Ma){M^ - M+){M^ - M_), (58) 

where we introduce the notation 

Q ^hi+ dawl, a = h^nox, c = X, (59) 

X= ^ — , M = , — — • (60) 

hi + daCw "ow hi + cj^da 

To study the modifications to the dynamics it is enough to expand each root in powers of the small parameter x 



Mo - ^x'\ M± ~ ± 
Pi 



1 + i (/3i - ai) f Ax + \ (5/3i - ai 



(61) 



where /3i = h^riQ and cii = g/(w^no). Since the parameter A and the momentum k appear only in the combination 
Afc, we have the symmetry property 

G^{uj,k) = G-^{uj,-k), (62) 

which will be useful in the final calculation of the Green functions G^{ijj,T{). In the radiation approximation, the 
integral in ([5^ produces 

G^f (^, R) = --^^ E \ - "^^'^ + *^«^Pfc"p) fc'^^e^"''' G\c^, A;), (63) 
(27r) ^ X 

where Ui = /r and R = |r — r' | . From now on we set i? = r in all denominators and understand that R — r — n-r' in 
the exponentials. We also neglect terms of order higher than 1/r. The cutoff O < Eqg defines the low energy domain 
of the effective theory. In this way we identify G^{uj, R) as 

G^(w,R) = 7-/ kdke'^^G^iuj,k). (64) 

The factor e^^^ forces us to close the integration contour in the upper half complex plane, choosing for example a 
semicircle with radius k — , picking up the poles in this region. Our description is valid only for momenta k <<Q.. 
According to Eqs. (pT|) . the pole at Mq corresponds to the momentum value |fco| = [Q^jT^I C^^- In the present 
approximation the contribution to the integral of this pole, together with the one of the semicircle in the upper half 
complex plane, can be neglected. The two remaining poles, which are the ones that contribute to the integral, are 
located at small displacements with respect to |fc-|-| — uquj << Eqg- In this way we take 

G\uj, k) = ^2^2QaMo (^A _ M+) {M^ - M_) ■ ^^^^ 



From the leading order expressions in (|6ip we conclude that the pole that contributes in the case A 
{lo + ie) uqAI^, while for A = — 1 it is — {u; + ie) noM_. The resulting integral is 



G^(w,R) = 



1 1 2nx 



-e 



iujnxR 



AnQ r n_ + n-|_ 

where we have considered that the dominant term in Mq yields gMq — 1. The refraction indices are 

nxiuj) = XnoMx. 
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-1 is 

(66) 
(67) 



The minus sign in ?i_ is because Af_ starts with a —1. Up to the order considered, the refraction indices are given 
by the expressions 



nx{uj) = no 



X 



(68) 



From Eq. ([55]) . and using Eqs. ([^ . (|27ll29p . ([50)1 and (pTj) . we can obtain the second order expansion for n\ 



nA(w) ~ 1 + A rii + (^^iS^ 



^2, 



with the real coefRcients ni and ^2 given by 

n-i = ("10 - ^lo) , = ^ [{aw - /3io) ("lo - 5/3io) + 4 (/3o2 - "02)] • 

Z o 

According to this the phase velocity is not 1. Due to the dispersive character of the background it becomes 



(69) 



(70) 



Vph{uj) ^ 1 - A (j_u?j rii + (^S_u?j (^{nif - 



Thus, the Green function in terms of space-time coordinates is 



G^(t,R) = — / cLu 



2n 



Q (n_ + n+) 



^ ^iujn\R^ — iujT 



duj 



Ji^ [l+A(|[^)ni + (|w)^n2] R^-iuT fjY^ 



where t = t — t' . 

If il ^ cx) the choice of the poles warrants the causal behavior of the Green function. But the frequency cutoff 
could introduce some violation of causality. To investigate this possibility, we compute the Fourier transform of the 
Green function to obtain its time dependent expression, by expanding the integrand in powers of ^ 



G'^(r,R)~— / cfw<j l + Ar^l(l + ^^^Ji?)wC- 
47^^ 



"20 - P02 — i {n2 + nf) Rlo + -n\R^ULi'^ 



^ <; 1 - anic (1 + RdR) du + ~e 



"20 - /3o2 - (?^2 + nl) ROr - -nlR^dl 



^2^sin|_^ (72) 



This shows that the main effect of the cutoff is to spread the propagating field around the light cone, within a 
wedge defined by ~ r ± 7r/2f2. Returning to G'^(a;,R), Eq. IpB]). we can characterize the effect of the cutoff in 



11 



the causal behavior of the Green function using the generahzed susceptibihty theorem [2l|, a generahzation of the 
Kramers-Kronig relations. Its real and imaginary parts as functions of the frequency lo are 



Re G^(uj,R) ~ 



1 + ( «20 - f3o2 - {n2 + nj) ROr - ^nlR^dl ) dl 



Im G^{lu,R) 



1 + ( "20 - Po2 - {n2 + nl) ROr - -nlR^dl ] dl 



cosLuR + XniS^ (1 + ROr) Or, sin {luR) j , 

(73) 

smojR — Xni£_ (1 + ROr) Or cosljR 

(74) 



To have a causal behavior they must satisfy the Kramers-Kronig relation 



1 



Im G{uj)\j^j^ ~ P \ duj 



,KeG{uj')-ReG{Vt) 



LO — LO 



(75) 



which gives 
ImG^c.)!^^. 



1 + r ( "20 - Po2 - [n2 + ni) ROr - -ntR'di ) 8^ 



/n 
duo 
-n 



,cos (io'R) — cos (flR) 



X no r. , ,sin(w'i?) -sin(f]i?)l 
~\nii{\ + RdR)dRP duo' ^ '- ^ '-) 



(76) 



For uo/Q, the integrals reduce to 



„ , , cos (lo'R) - cos (nR) lo 
PI duj' ^ ^ ^ ^ ~2(l-coswi?)-cos(17i?) 



n lo' - LO 

n 



n + [{flR) (cosfii?) - sinflR] 



,sin (lo'R) - sin^QR) 



P I duo' ^ "'^^ — "'"^"""^ ~ 2 (1 - sinwi?) ^ cos (flR) 



TT + [{nR) {cosflR) - sinnR] 



sin wi?, 
cos loR. 



(77) 



Furthermore, in the case of a radiation field ilR ^ 1 and hence the factors cosili? and sinili? become strongly 
oscillating, nullifying the contributions of the terms where they appear (which also have a factor (tj/ri)" ,with n> 1). 
Thus we can take 



, cos{uo'R)-cos{nR) , sin(^'j?)-sin(»j;) 
P I dio ; ~ —-KsmuoR, PI duj ; ~ tt cos wit. 

-n Lo' -uo Lo' -bj 



(78) 



Replacing these integrals in ([75]) . we finally get that if Re G [uo) is given by Eq. ([75]) . the imaginary part of the Green 
function, Im G^(cij), must be 



Im G\io)\^^ = — 



1 

47rr 



1 + ( "20 - /3o2 - {n2 + nl) ROr - -nfR^dl ) dl 



sin ujR — Ani^ (1 + ROr) Br cos ujR J , 

(79) 



which coincides with Eq. (|74p . obtained by direct computation. 



V. FINAL REMARKS 



To summarize, in the preceding sections we have presented a general description for a large class of effective models 
for the electromagnetic field incorporating dynamical corrections motivated by QG and leading to departures from 
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standard physics. The main features characterizing the models to which such a description is apphcable are: (1) 
the vahdity of gauge invariance and charge conservation, (2) the use of standard commuting space-time coordinates 
together with the corresponding Fourier transform methods (which is not the case of Double Special Relativity models, 
for example), (3) the assumption that effective field theories constitute an appropriate tool for describing the low energy 
behavior of remnant effects which could arise from quantum gravity, (4) the assumption that low energy dynamics 
is linear in the potential field, (5) the inclusion of non-local effects via the operator character of the generalized 
susceptibilities. This description makes it also possible to include anisotropic corrections in the constitutive relations, 
for example via additional non-dynamical tensors arising from spontaneous Lorentz symmetry breaking, a case which 
is not considered in this work. 

The proposed formalism is quite similar to the usual electrodynamics in a medium, except for the non-local character 
of the effective QG corrections, mirroring the granularity of the space-time induced by quantum fluctuations of the 
metric. This feature leads to an electrodynamics with non-local constitutive relations, which contains terms connecting 
D and H with both E and B. Thus the QG modifications can be modelled by a dispersive bianisotropic medium, 
where the propagation of the electromagnetic field is characterized by a refraction index, whose first order term in 
the perturbative parameter ^ is directly related to vacuum birrefringence. We have considered the models from the 
point of view of active transformations, i.e. observable Lorentz symmetry violations associated with boosts in a given 
reference frame 

The effective models correspond in fact to high order theories. Hence we used an adequate generalization of the 
Noether theorem to find the energy-momentum tensor to second order in the LIV parameter ^. Next we determined 
the density of energy and momentum carried by the electromagnetic field, for which we give explicit expressions to 
order ^. They acquire a simple form when written in terms of the fields E, B and D, H , which shows the convenience 
of the latter for describing the dynamics, in an analogous way to the usual electrodynamics in media. 

This theory is valid only for low energies. We have also studied the consequences of this fact by using an explicit 
cutoff il ^ Eqc ^ C^^- III i'ACt, the results in Eqs. ([77]) and ([75]) show that the introduction of the cutoff does not 
produce any significant causality violation in the radiation regime {flR >>> 1) because the expected modifications 
proportional to (w/il)" in (|77p . which are the subject of possible signals of new physics in these approaches, are 
further suppressed by highly oscillating terms proportional to sin(r2i?), cos(f2i?) thus nullifying the impact of causality 
violation upon the corresponding observational effects. The most outstanding manifestation of the cutoff is a spreading 
of the propagation of the electromagnetic field around the light cone. In fact there are two sources for such a spreading. 
One is due to the cutoff and the other arises from the dispersive character of the effective medium, which leads to an 
cj-dependent phase velocity. The relation between both effects depends on the relative value of fl and In any 
case, for distances large enough from the source {ujR ^1), the dispersive effect will finally dominate. 

There remains to discuss the causal behavior of the full theory. There are two possible sources of acausality. One 
is related to the dispersive character of the effective medium, while the other is related to the existence of velocities 
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V > 1 that leads to photons propagating to the past in highly boosted reference frames, and hence to the possibility 
of acausal loops. This issue is beyond the scope of the present work, and will be discussed in detail elsewhere. 

Acknowledgements 

R.M. acknowledges partial support from CONICET- Argentina. L.F.U is partially supported by projects 
CONACYT-40745F, CONACYT-47211F and DGAPA-UNAM-IN104503-3. 



[1] Amelino-Camelia, G., Ellis, J., Mavromatos, N. E., Nanopoulos, D. V. and Sarkar, S., Nature (London) 393 (1998) 763. 
[2] Gambini, R. and PuUin, J., Phys. Rev. D59 (1999) 124021; Gambini, R. and Pullin, J., in Proceedings of the Second 

Meeting on CPT and Lorentz Symmetry, ed. Kostelecky, V.A. (World Scientific, 2002). 
[3] Alfaro, J., Morales- Tecotl, H.A., and Urrutia, L. F., Phys. Rev. Letts. 84 (2000) 2318; Alfaro, J., Morales-Tecotl, H.A., 

and Urrutia, L. F., Phys. Rev. D65 (2002) 103509; Alfaro, J., Morales-Tecotl, H.A., and Urrutia, L. P., Phys. Rev. D66 

(2002) 124006. 

[4] Biller, S.D., et al, Phys. Rev. Lett. 83 (1999) 2108; Bertolami, O. and Carvalho, C. S., Phys. Rev. D61 (2000) 103002; 
Lammerzhal C. and Borde C., in Gyros, Clocks, Interferometers,... : Testing Relativistic Gravity in Space, Lecture Notes 
in Physics 562 (Springer, Heidelberg 2001); Haugan , M. P. and Lammerzhal, C., in Gyros, Clocks, Interferometers,... 
: Testing Relativistic Gravity in Space, Lecture Notes in Physics 562 ( Springer, Heidelberg 2001); Gleiser, R. J., and 
Kozameh, C. N., Phys. Rev. D64 (2001) 083007; Brustein, R., Eichler, D., and FoflFa, S., Phys. Rev. D65 (2002) 105006; 
Konopka, T. J. and Major, S. A., New J. Phys. 4 (2002) 57; Sudarsky, D., Urrutia, L. F. and Vucetich, H., Phys. Rev. 
Lett. 89 (2002) 231301; Alfaro, J. and Palma, G., Phys. Rev. D67 (2003) 083003; Lambiase, G. and Singh, P., Phys. Letts. 
B565 (2003) 27; Stecker, F. W., J. Phys. G29 (2003) R47; Ellis, J., Mavromatos, N.E., Nanopoulos, D.V., and Sakharov, 
A.S., Astron. Astrophys. 402 (2003) 409; Ellis, J., Mavromatos, D.V., and Sakharov, A.S., Astropart. Phys. 20/6 (2004) 
669; Ellis, J., Mavromatos, N.E., Nanopoulos, D.V., and Sakharov, A.S., Nature (London) 428 (2004) 386. 

[5] Jacobson, T., Liberati, S. and Mattingly, D., Nature (London) 424 (2003) 1019; ibid., Phys. Rev. D67 (2003) 124011; 
ibid., Phys. Rev. D66 (2002) 081302; Jacobson, T., Liberati, S., Mattingly, D. and Stecker, F. W., Phys. Rev. Letts. 93 
(2004) 021101. 

[6] CoUaday, D., and Kostelecky, V.A., Phys. Rev. D55 (1997) 6760; Phys. Rev. D58 (1998) 116002. 

[7] For a review see for example Proceedings of the Meeting on CPT and Lorentz Symmetry, ed. Kostelecky, V.A., (World 
Scientific, 1999); Will, C. M., The Confrontation between General Relativity and Experiment (Living Reviews in Relativity 

4, 2001); Proceedings of the Second Meeting on CPT and Lorentz Symmetry, ed. Kostelecky, V.A. (World Scientific, 2002). 
[8] Thiemann, T., Class. Quant. Grav. 18 (2001) 2025; Thiemann, T., and Winkler, O., Class. Quant. Grav. 18 (2001) 2561; 

ibid.. Class. Quant. Grav. 18 (2001) 4629; ibid.. Class. Quant. Grav. 18 (2001) 4997; Sahlmann, V., Thiemann, T., and 
Winkler, O., Nucl. Phys. B606 (2001) 401. 
[9] Ellis, J., Mavromatos, N.E., and Nanopoulos, D.V., Gen. Rel. Grav. 32 (2000) 127; Phys. Rev. D 61 (2000) 027503; 
Ellis, J., Farakos, K., Mavromatos, N.E., Mitsou, V.A. and Nanopoulos, D.V., Astrophysical Jour. 535 (2000) 139; Ellis, 
J., Mavromatos, N.E., Nanopoulos, D.V., and Volkov, G., Gen. Rel. Grav. 32 (2000) 1777. 
[10] Myers, R.C. and Pospelov, M., Phys. Rev. Letts. 90 (2003) 211601 ; Perez, A. and Sudarsky D., Phys. Rev. Letts. 91 

(2003) 179101. 

[11] R. Montemayor and L. F. Urrutia, Phys. Rev. D72, 045018 (2005); ibid., AIP Conf. Proc. 758, 81 (2005); ibid., Phys. 
Lett. B606, 86 (2005). 

[12] Collins, J., Perez, A., Sudarsky, D., Urrutia, L. F. and Vucetich, H., Phys. Rev. Lett. 93 (2004) 191301 . 
[13] Alfaro, J., Phys. Rev. Lett. 94 (2005) 221302 ; ibid. Phys. Rev. D72 (2005) 024027. 
[14] Kostelecky, V. A. and Lehnert, R., Phys. Rev. D63 (2001) 065008. 

[15] See for example Rovelh, C. and Speziale, S., Phys. Rev. D67 (2003) 064019; Dowker, P., Henson J. and Sorkin, R. D., 
Mod. Phys. Lett. 19A (2004) 1829; Magueijo, J. and Smolin, L., Phys. Rev. Lett. 88 (2002) 190403; Magueijo, J. and 
Smolin, L., Phys. Rev. D67 (2003) 044017; Girelli, F. and Livine, E., Physics of deformed special relativity: relativity 
principle revisited, |gr-qc/0412004, Chryssomalakos, C. a nd Ok on, E., Int. J. Mod. Phys 13D (2004) 2003; Smolin, L., 
Falsifiable predictions from semiclassical guantum gravity, |hep-th/0501091 Corichi, A., Sudarsky, D., Int. J. Mod. Phys. 
D14 (2005) 1685. 

[16] Amelino-Cameha, G., Phys. Lett. B510 (2001) 255; ibid., Int. J. Mod. Phys. Dll, (2002) 35. 

[17] For a review see for example: Amelino-Camelia, G., New J. Phys. 6 (2004) 188 , Vucetich, H., Testing Lorentz invariance 
violation in quantum gravity theories, gr-qc/0502093 ; Mattingly, D., Living Rev. Rel. 8 (2005) 5; Jacobson, T., Liberati, 

5. and Mattingly, D., Annals Phys. 321 (2006) 150. 

[18] Landau, L.D. and Lifshitz, E. M., The Classical Theory of Fields, Third Edition, Addisson- Wesley, Chapter 10. 



14 



[19] Gonzalez, B., Martinez, S. A., Montemayor, R. and Urrutia, L. F., Journal of Physics: Conference Proceedings Series 24 

(2005) 58. 

[20] Martinez, S. A., Montemayor, R. and Urrutia, L. F., Phys. Rev. D74 (2006) 065020. 
[21] Landau, L.D. and Lifshitz, E. M., Physique Statistique, editions MIR, Moscou 1984. 



